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We consider a magnetized degenerate gas of fermions as the matter source of a homo- 
geneous but anisotropic Bianchi I spacetime with a Kasner-like metric. We examine the 
dynamics of this system by means of a qualitative and numerical study of Einstein- 
Maxwell field equations which reduce to a non-linear autonomous system. For all initial 
conditions and combinations of free parameters the gas evolves from an initial anisotropic 
singularity into an asymptotic state that is completely determined by a stable attractor. 
Depending on the initial conditions the anisotropic singularity is of the "cigar" or "plate" 
types. 
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1. Introduction 

Intense magnetic fields may play an important role in astrophysical objects, such 
as compact object like neutron or white dwarf stars. Weaker magnetic fields ^ 
may also be important in connection with the magnetic plasma effects associated 
with local anisotropy in self-gravitating and collision-less systems (e.g., galactic ha- 
los of fermionic dark mater l™J and net radiation flow entering (or leaving) the gas 
clouds. Strong magnetic fields can also arise as plausible explanation behind the ori- 
gin of anisotropy in highly dense systems, like solid core, exotic phase transitions, 
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pion condensation or anisotropic compact object by ultra-strong self-magnetization 
(B > 10 14 ) Gauss. The latter phenomenon motivates the present work, as a the- 
oretical explanation behind th e aniso tropy of stresses in these systems. This is 
connected with previous 

wor k | 4 | 5 | 7 | 6 | in which 

we studied the thermodynamical 
properties of gases of electrons and vectorial bosons using the Electroweak model. 
We showed in that work that such a system was characterized by an anisotropic 
energy-momentum tensor in which the pressure transverse to the magnetic field 
(p_i_) was (in general) different from that parallel to the field (py ^ p±). Under a 
Newtonian framework, we found that the vanishing of p± would drive the system 
to gravitational collapse. We also discussed if critical values of the magnetic field 
exist, allowing for the vanishing of this pressure and its subsequent collapse. By 
assuming that the model could describe a Newtonian compact object, we concluded 
that the final state of such a star would never be a magnetar of the type of Duncan 
and Tompson^. Instead, such a state could be a strange star or a "black cigar-like" 
structure. 

However, the proper study of highly magnetized sources of possible compact 
objects requires that we incorporate the non-linear effects of Einstein's General 
Relativity theory, allowing us to examine the effects of the gas magnetization in 
the collapse of a magnetized self-gravitating Fermi gas. Ideally, we should consider 
minimally realistic rotating configurations, but this would require elaborate hydro- 
dynamical numerical techniques. Instead and as a first step, we consider simple 
and mathematically more tractable, even if extremely idealized, configurations that 
could reveal (at least) some of the properties of the dynamical evolution of magne- 
tized gases in General Relativity. We undertake such a study in the present article 
by considering what could be the simplest non-static geometry compatible with 
the anisotropies of a magnetized gas: the homogeneous but anisotropic Bianchi I 
spacetime in a Kasner-like coordinate basis. 

2. Magnetized Fermi gas for a Kasner metric. 

The Bianchi I spacetime can be described by the Kasner-like metric element EES 

ds 2 = -c 2 dt 2 + A 2 {t)dx 2 + B 2 {t)dy 2 + C 2 {t)dz 2 , (1) 

corresponding to an homogeneous, but anisotropic, spacetime with zero spacial cur- 
vature and having all relevant quantities depending only on time. We consider as a 
source of ([1]) a degenerated magnetized Fermi gas whose energy-momentum tensor 
is that obtained by A. Perez Martinez and H.Perez Rojas^. 

The energy-momentum tensor for this gas in a comoving frame with coordinates 
x = ct, x, y, z and 4- velocity u a = 5q is given by 

T b a = (U + P)u a u b + P5- + U a b , P = p-^BM, (2) 

where U is the total matter-energy density, B the magnetic field and M. the magneti- 
zation of the gas, P is the isotropic pressure and is the traceless anisotropic pres- 
sure tensor, which for flU) takes the form U b = diagTl, II, —211, 0] with n = ~BM. 
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Thus, p — T* is the stress tensor component (pressure) in the direction of the 
magnetic field. 

The full equation of state for this gas takes the following form 

p= \pT p {p,n), BM = \pT M (p,n), U = \pT ff (p,n). (3) 

where the T functions are given further ahead, /3, the magnetic field normalized 
with a suitable critical value B c , and A are given by: 

B inc 2 

If we consider our gas to be made of by electrons, then m is the electron mass, m e , 
while A c = h/mc is the Compton wavelength and B c = m 2 c 3 /eH is the so-called 
critical magnetic field 0. In this limit of strong magnetic field, the functions T(/3, fj,) 
take the form: 



^ a n{^"n b n C n ), Yp ^ Q? n C n , Yp -\- 2 ^ Oi n (X n (5) 

n— n— n— 



where a„ = ^vV 2 - 1 - 2n/3, 6„ = ln[(// + a„//i)/ % /l + 2n/3J, c„ = 2n[3b„, a n = 2 — 
Son, and /i = /i e /(m e c 2 ) chemical potential of the fermions, s = I[^(B C /B)(fi 2 — 1)] 
corresponds to the maximum Landau level for a given Fermi energy and magnetic 
field strength, while I[X] denotes the integer part of its argument X. 

Although {1} is obviously inappropriate to examine the magnetized g as as a 
source of a compact object, we will consider in our analysis the conditions pre- 
vailing inside a very high density compact object in which we can ignore the gas 
temperature. 



3. Einstein Maxwell field equations. 

The comoving 4- velocity u a — Sq in fl} yields the expansion scalar 8 = u a - a and the 
traceless shear tensor a a b — ^(a-,b) ~ (l/3)0/i a 6- In order to obtain a self-consistent 
system of ordinary differential equations suitable for numerical integration, it is con- 
venient to eliminate second order derivatives of the metric functions by expressing 
the coupled Einstein-Maxwell field equations G£ = kT^, F ab $ = and F[ a j,- C ] = 0, 
as well as the balance equations T ab -j, = 0, in terms of 8 and cr£. Following this 
approach and after tedious algebraic manipulations, we can reduce the field and 
balance equations to the following self-consistent system of 4 independent ordinary 



a (e.g. for the electron B% = 4, 41 X 10 13 Gauss) 
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differential equations: 



(U+p--BM)e-BME 



U 

£ = ~r:BM -HI 

e 



{bm + -(u - P ))-e 2 

/3(3S - 26) 



(6a) 
(6b) 
(6c) 
(6d) 



where £ = erf and a dot denotes derivative with respect to ct. From ([5]), the state 
variables C/, p and S.M depend on the normalized magnetic field and chemical po- 
tential, ft and fj,. 

It is convenient to assume the zero order approximation (n = 0), so that we 
would be considering all electrons in the basic Landau level Under this approxi- 
mation the expansions of the T(f3,n) = r(/i) functions in (|5|) simplify considerably: 



(ED- 



ao — bo and T u — clq + bo, where an, &o are given by setting n = in 

Introducing the following new variables and dimensionless functions: 

kA/3 



„ _ e _d_ s 

3' Bcdt' ff' 



n = 



3ff 2 



« = ( 7 ) 

-no 



where Hq is a constant with inverse length units. Einstein-Maxwell equations ([6]) 
become the following dimensionless system 



ft' = 
S' = 

p> = 

supplemented by 

w = 



2 { 1 S 



1 + 2 In 



(2 - S) ny/tj? - 1 - 2(1 + S) In L + vV - l) 



n, 



(2 - 5) In (ji + ^W^l) - SSfi^i^T 



2fi 2 



V7 



{-3 + [/i^/i 2 - 1 + 2 In ffi + vV - l)] ft} 



Ti. 



(8a) 
(8b) 

(8c) 
(9) 



where a prime denotes derivative with respect to the dimensionless time r and we 
have chose a length scale characterized by setting (3Hq/kX) = 1, so that (from 
0), we obtain 



ft 



nn 2 . 



(10) 



Since we have k\ = 0.749 x 10 24 cm 2 for electrons, this choice leads to 1/Hq = 
2 x 10 12 cm, which is much smaller than the cosmological Hubble radius. E| 



b It is of the order of the distance from the Earth to the Sun. 
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4. Qualitative and numerical analysis. 

The system ([8]) defines a three dimensional phase space (£7, .S', /x) with £2 > and 
fi > 1. The planes £7 = and /i = 1 are invariant subsets corresponding to a Kasner 
vaccum. Besides these subsets, we also have: 

(11a) 

2b (-l + 2/i 2 + 2a ) 1 

2& )ao + (-l + 2M 2 )&o-6M 2 (l- M 2 )' J ' 

(lib) 

III = { A t = 1.42, S = 0.34, £2 = 0.42} (11c) 

where ao, &o are given by setting n = in ([5]). The curves I and II are saddles, while 
III is a stable attractor. As shown in figure 1, some trajectories start in a singular 
state characterized by £2 — > oo while others start in a singular Kasner vacuum 
£2 = 0, though all terminate in the attractor III. Depending on initial conditions, 
some trajectories approach eras of Kasner vacum I, II. The asymptotic state is given 
by III, which (from (fT0|) ) corresponds to (3 — > 0.42 Ti 2 . Numerical examination of ((8]) 
and ([9]) reveals that (3 has a an asymptotic power law decay in terms of the physical 
time t. Trajectories that initiate in a singular Kasner vacum denote an unphysical 
evolution passing from £2 = to the value £2 = 0.42 given by the attractor III. The 
most physically motivated trajectories correspond to those in which £2 decreases 
monotonically to the asymptotic value in III, hence the normalized magnetic field 
(3 also exhibits monotonous decay approximately following a power law scaling law 
(in terms of the physical time t) . 

Since the physically interesting range of magnetic field intensities in a magne- 
tized electron gas is B > B c (or (3 > 1), it is important to examine the behavior of 
the trajectories near the initial singular states. The nature of the initial singularity 
can be appreciated by looking at the behavior of the scale factors A, B. C in ^ 
in the limit t — > — oo. The initial singularity (whether a Kasner vaccum or not) is 
generically anisotropic and can be of the "cigar" or "plate" types. In the former 
case we have C — > oo with diverging £2, while both A, B tend to zero, forming a 
singular line along the z direction parallel to the magnetic field. 



5. Conclusion 

We have examined the evolution of magnetized degenerate gas of fermions, consid- 
ering the case in which the fermions are electrons in their basic Landau level (n = 
in ©). Einstein-Maxwell equations for a Bianchi I spacetime yield the dynamical 
system ([8]) . Initial conditions exist for an evolution starting in a line singularity 
with £2 (and (3) decaying monotonously to the asymptotic value in the attractor III 
This type of evolution is similar to that emerging in previous Newtonian studies 
However, we have also the possibility of a non-vaccum initial singular state that is 
not unphysical (diverging £2) and is characterized by a "plate" type of singularity, 
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Fig. 1. Phase space (Q, S, (j,). Notice how all trajectories evolve towards the attractor III. Some 
curves approach states of Kasner vacum (the saddles Q = 0, I and II). Some curves have an 
initial singular vacum state (III) and some initiate at f! —* oo. Singular state are characterized by 
anisotropic singularities. 



which is a countcr-intutive type of anisotropy that lacks a Newtonian equivalent. A 
more detailed and extensive study of the singular states in this gas, as well as a less 
restrictive form of the equation of state and the case of a gas of neutrons (instead 
of electrons), is currently under investigation. 
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